Dimensionality and dynamics in the behavior of C. elegans 
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A major challenge in analyzing animal behavior is to discover some underlying simplicity in 
complex motor actions. Here we show that the space of shapes adopted by the nematode C. elegans 
is surprisingly low dimensional, with just four dimensions accounting for 95% of the shape variance, 
and we partially reconstruct 'equations of motion' for the dynamics in this space. These dynamics 
have multiple attractors, and we find that the worm visits these in a rapid and almost completely 
deterministic response to weak thermal stimuli. Stimulus-dependent correlations among the different 
modes suggest that one can generate more reliable behaviors by synchronizing stimuli to the state of 
the worm in shape space. We confirm this prediction, effectively "steering" the worm in real time. 



INTRODUCTION 

The study of animal behavior is rooted in two diver- 
gent traditions. One approach creates well-controlled 
situations, in which animals are forced to choose among 
a small discrete set of behaviors, as in psychophysical 
experiments [ij. The other, taken by ethologists 2], de- 
scribes the richness of the behaviors seen in more natural 
contexts. One might hope that simpler organisms pro- 
vide model systems in which the tension between these 
approaches can be resolved, leading to a fully quantita- 
tive description of complex, naturalistic behavior. 

Here we explore the motor behavior of the nematode, 
Caenorhabditis elegans, moving freely on an agar plate 
0, 0i S, H 0] ■ Though lacking the full richness of a nat- 
ural environment, this unconstrained motion allows for 
complex patterns of spontaneous motor behaviors Q, 
which are modulated in res pon se to chemical, thermal 
and mechanical stimuli 0, llOl [llj . Using video mi- 
croscopy of the worm's movements, we find a low di- 
mensional but essentially complete description of the 
macroscopic motor behavior. Within this low dimen- 
sional space we reconstruct equations of motion which 
reveal multiple attractors — candidates for a rigorous def- 
inition of behavioral states. We show that these states 
are visited as part of a surprisingly reproducible response 
of C. elegans to small temperature changes. Correlations 
among fluctuations along the different behavioral dimen- 
sions suggest that some of the randomness in the behav- 
ioral responses could be removed if sensory stimuli are 
delivered only when the worm is at a well defined initial 
state. We present experimental evidence in favor of this 
idea, showing that worms can be "steered" in real time 
by appropriately synchronized stimuli. 



EIGENWORMS 

We use tracking microscopy with high spatial and tem- 
poral resolution to extract the two-dimensional shape 



of individual C. elegans from images of freely moving 
worms over long periods of time (Fig la; see Methods). 
Variations in the thickness of the worm are small, so we 
describe the shape by a curve that passes through the 
center of the body (Fig lb). We measure position along 
this curve (arc length) by the variable s, normalized so 
that s = is the head and s = 1 is the tail. The posi- 
tion of the body element at s is denoted by x(s), and we 
sample this function at = 100 equally spaced points 




FIG. 1: Describing the shapes of worms, (a) Raw image in 
the tracking microscope, (b) We process the image in (a) and 
pass a curve through the center of the body. The black circle 
marks the head, (c) We measure distance along the curve 
(arclength s) in normalized units, and at each point we define 
the tangent t(s) and and normal n(s) to this curve. The 
tangent points in a direction 0{s), and variations in this angle 
correspond to the curvature k(s) = d9{s)/ds. (d) We rotate 
all images so that the average angle 9 is zero; therefore d{s) 
provides a description of the worms shape that is independent 
of our coordinate system, and intrinsic to the worm itself. 
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along the body. These variables provide an essentially 
complete description of the motor output. 

We analyze the worm's shapes in a way intrinsic to its 
own behavior, not to our arbitrary choice of coordinates. 
The intrinsic geometry of a curve in the plane is defined 
by the Frenet equations l3, 



(a) 



dx{s) 

ds 
di{s) 

ds 



K{s)n{s), 
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where t(s) is the unit tangent vector to the curve, n(s) is 
the unit normal to the curve, and k{s) is the scalar cur- 
vature (Fig lb). If the tangent vector points in a direc- 
tion 0{s), then k{s) = dO{s)/ds. Curvature as a function 
of arc length, k(s), thus provides a "worm-centered" de- 
scription, but in practice this involves taking two deriva- 
tives and thus is noisy. As an alternative, we describe 
the curve by 9{s), but remove the dependence on our 
choice of coordinates by rotating each image so that the 
mean value of along the body always is zero (Fig Ic); 
this rotated version of 6{s) contains exactly the same 
information as k{s). 

Although the worm has no discrete joints, we expect 
that the combination of elasticity in the worm's body 
wall and a limited number of muscles will lead to a lim- 
ited effective dimensionality of the shape and motion. 
In the simplest case, the relevant low dimensional space 
will be a Euclidean projection of the original high dimen- 
sional space. If this is true, then the covariance matrix 
of angles, C(s,s') = {{0{s) - {0)){O{s') - {9))) will have 
only a small number of nonzero eigenvalues. Figure 2a 
shows the covariance matrix, and its smooth structure 
is a strong hint that there will be only a small number 
of significant eigenvalues; this is shown explicitly in Fig 
2b. Quantitatively, if we measure the total variance in 
angle along the body, then over 95% of this variance can 
be accounted for by just four eigenvalues. 

Associated with each of the eigenvalues is an eigen- 
vector u^{s), sometimes referred to as a 'principal com- 
ponent' of the function 9{s). If only K — A eigenvalues 
are significant, then we can write the shape of the worm 
as a superposition of 'eigenworm' shapes. 
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FIG. 2: Covariance of shape fluctuations and eigenworms. 
(a) The covariance matrix of fluctuations in angle C(s, s') is 
computed from 60, 000 images across eight diflterent worms. 
The inhomogeneity along the diagonal shows that the nor- 
mal modes of the motion are not sinusoidal but the smooth 
structure of C'{s,s') means that a small number of modes 
are significant, (b) We find the eigenvalues of C(s,s') and 
compute a]f, the fraction of the total variance captured by 
keeping just K modes (see Methods); 95% of the variance is 
accounted for by the first four modes, (c) Associated with 
each of the dominant modes is an eigenvector describing an- 
gle vs. position along the worm; we refer to these as the 
eigenworms itp(s). The population-mean eigenworms (red) 
are highly reproducible across individual worms (black), (d) 
In response to strong thermal stimuli the worm's shape re- 
mains within the four dimensional subspace identified from 
analysis of spontaneous crawling. Worm images are recorded 
at times synchronized to a heat pulse (see Methods) and we 
display cr^ aUgned with the stimulus. Though the shapes 
vary systematically in response to these stimuli, the accu- 
racy of our reconstructions does not and using all four modes 
we continue to account for 95% of the shape variance. 



9{s) 



K 



(3) 



where the four variables {a^} are the amplitudes of 
motion along the different principal components, = 
X^s We see in Fig 2c that these modes are 
highly reproducible from individual to individual. 

Thus far we have considered only worms moving in 
the absence of deliberate sensory stimuli. Do the worms 
continue to move in just a four dimensional shape space 
when they respond to strong inputs? To test this, we 



delivered intense pulses of heat (see Methods) , which are 
known to trigger escape responses Q . We see in Fig 2d 
that we still account for ^ 95% of the shape variance 
using just four modes, even though the distribution of 
shapes during the thermal response is very different from 
that seen in spontaneous crawling. We conclude that our 
four eigenworms provide an effective, low dimensional 
coordinate system within which to describe C. elegans 
motor behavior. 
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FIG. 3: Motions along the first two eigenworms. (a) The 
joint probabihty density of the first two amphtudes, p{ai, 02), 
with units such that (ai) = (02) = 1. The ring structure sug- 
gests that these modes form an oscillator with approximately 
fixed amplitude and varying phase = tan^^ (— a2/ai). (b) 
Images of worms with different values of show that variation 
in phase corresponds to propagating a wave of bending along 
the worms body, (c) Dynamics of the phase (j){t) shows long 
periods of linear growth, corresponding to a steady rotation 
in the {01,02} plane, with occasional, abrupt reversals, (d) 
The joint density p{\v\, \u!\). The phase velocity to = d(j>/dt 
in shape space predicts the speed at which the worm crawls 
on the agar. The crawling speed was defined as the time 
derivative of the worm's center of mass. 



WHAT DO THE MODES MEAN? 

The projection of worm shapes onto the low- 
dimensional space of eigenworms provides a new and 
quantitative foundation for the classical, qualitative de- 
scriptions of C. elegans behavior Q . The first two modes 
are sinuous (although not exactly sinusoidal) oscillations 
of the body shape (Fig 2c); they form a quadrature 
pair, so that different mixtures of the two modes cor- 
respond to different phases of a wave along the body. 
Indeed, the probability distribution of the mode ampli- 
tudes, p(ai,a2), shows a ring of nearly constant ampli- 
tude (Fig 3a) . Sampling images around this ring reveals 
a traveling wave along the body (Fig 3b). There are 



relatively long periods of time where the shape changes 
by a continuous accumulation of the phase angle (Fig 
3c), and the speed of this rotation predicts the speed at 
which the worm crawls (Fig 3d). 

In contrast to the first two modes, the third 
mode M3(s) contributes to a nearly constant curvature 
throughout the middle half of the body (Fig 2c). The 
distribution of the mode amplitude 03 has a long tail 
(Fig 4a), and body shapes chosen from these tails (Fig 
4b) exhibit the O configuration classically identified with 
turning behavior [3] . Large amplitudes of 03 also corre- 
spond to regions of high curvature in the worm trajectory 
along the agar (Fig 4c). 

The connections between mode amplitudes and the 
motion of the worm along the agar — as in Figs and 
|3h — are genuine tests of the functional meaning of our 
low dimensional description. Quite explicitly, our anal- 
ysis of worm shapes is independent of the extrinsic co- 
ordinates and hence our definition of modes and ampli- 
tudes is blind to the actual position and orientation of 
the worm. Of course, in order to move the worm must 
change shape, and our description of the shape in terms 
of mode amplitudes captures this connection to move- 




(c) 



FIG. 4: Motions along the third eigenworm. (a) The dis- 
tribution of amplitudes ^(03), shown on a logarithmic scale. 
Units are such that (ai) = 1, and for comparison we show the 
Gaussian distribution; note the longer tails in p{a3). (b) Im- 
ages of worms with values of 03 in the negative tail (left), the 
middle (center) and positive tail (right). Large amplitudes of 
03 correspond to Sl-like shapes, (c) A two minute trajectory 
of the center of mass sampled at 4 Hz. Periods where |a3 1 > 1 
are colored red, illustrating the association between turning 
and large displacements along this mode. 



4 



ment. Thus, to crawl smoothly forward or backward the 
worm changes its shape by rotating clockwise or counter- 
clockwise in the plane formed by the mode amplitudes 
oi and 02; the speed of crawling is set by the speed of 
the rotation. Similarly, to change direction the worm 
changes shape toward larger magnitudes of the mode 
amplitude 03, and we see this connection even without 
defining discrete turning events. 



ATTRACTORS AND BEHAVIORAL STATES 

The eigenworms provide a coordinate system for the 
postures adopted by C. elegans as it moves; to describe 
the dynamics of movement we need to find equations of 
motion in this low dimensional space. We start by focus- 
ing on the plane formed by the first two mode amplitudes 
oi and 02. Figure [3] suggests that within this plane the 
system stays at nearly constant values of the radius, so 
that the relevant dynamics involves just the phase angle 
(j){t). To account for unobserved and random influences 
these equations need to be stochastic, and to support 
both forward and backward motion they need to form a 
system of at least second order. Such a system of equa- 
tions would be analogous to the description of Brownian 
motion using the Langevin equation ^ISt . il6i | . Thus we 
search for equations of the form 



d4>{t) 

dt 
duj{t) 

dt 



= Fm),u;it)]+am),u{tMt). (4) 



Here F[(l){t),uj{t)] defines the average acceleration as a 
function of the phase and phase velocity, by analogy to 
the force on a Brownian particle. The noise is charac- 
terized by a random function 77 (i) which we hope will 
have a short correlation time, and we allow the strength 
of the noise a[<j){t),uj{t)] to depend on the state of the 
system, by analogy to a temperature that depends on 
the position of the Brownian particle. 

In Figure [5|i we show our best estimate of the mean 
acceleration F[(l),uj] (see Methods for details). Once we 
know F, we subtract this mean acceleration from the 
instantaneous acceleration to recover trajectories of the 
noise, and the correlation function of this noise is shown 
in Fig 5b. The correlation time of the noise is short, 
which means that we have successfully separated the dy- 
namics into two parts: a deterministic part, described by 
the function _F[0,a;], which captures the average motion 
in the {ai, 02} plane and hence the relatively long peri- 
ods of constant oscillation, and a rapidly fluctuating part 
r]{t) that describes "jittering" around this simple oscil- 
lation as well as the random forces that lead to jumps 
from one mode of motion to another. 

We can imagine a hypothetical worm which has the 
same deterministic dynamics as we have found for real 



-<tl(t)ll(t+T)> 
-<0j(t)u3(t+x)> 




FIG. 5: Reconstructing the phase dynamics, (a) The mean 
acceleration of the phase as a function of the instantaneous 
phase and phase velocity, F{uj, (j)) in Eq (5). (b) The correla- 
tion function of the noise in the phase dynamics, {r]{t)rj{t+T)) 
with ri{t) defined from Eq (5). The correlations are confined 
to very short times, especially when compared with the cor- 
relations of the phase velocity itself, (c) Trajectories in the 
deterministic dynamics. A selection of early-time trajectories 
are shown in black. At late times these same trajectories col- 
lapse to one of four attractors (red): forward and backward 
crawling and two pause states, (d) Attracting basins. Ini- 
tial conditions covering the {(p, u)} plane were evolved to long 
times and each initial phase point was colored corresponding 
to its asymptotic attractor. Note that the basins are asym- 
metric so that some forward motions fall into the reversal or 
pause attractor. This highlights the difference between a dy- 
namical definition of behavioral state and phenomenological 
one based on an instantaneous value such as sgn{uj{t)). 



worms, but no noise. We can start such a noiseless worm 
at any combination of phase and phase velocity, and fol- 
low the dynamics predicted by Eqs but with a = 0. 
The results are shown in Fig[St. The dynamics are di- 
verse on short time scales, depending in detail on the 
initial conditions, but eventually all initial conditions 
lead to one of a small number of possibilities (Fig[5jl): 
either the phase velocity is always positive, always neg- 
ative, or decays to zero as the system pauses at one of 
two stationary phases. Thus, underneath the continu- 
ous, stochastic dynamics we find four discrete attractors 
which correspond to well defined classes of behavior. 



PAUSE STATES AND REPRODUCIBILITY 

The behavior of C. elegans, particularly in response 
to sensory stimuli, traditionally has been characterized 
in probabilistic terms: worms respond by changing the 
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probability of turning or reversing [B, H This 
randomness could reflect an active strategy on the part 
of the organism, or it could reflect the inability of the 
nervous system to distinguish reliably between genuine 
sensory inputs and the inevitable background of noise. 
Our ability to describe motor behavior more completely, 
with high time resolution, offers us the opportunity to 
revisit the "psychophysics" of C. elegans. 

We consider the response to brief (75 ms), small 
(AT ~ 0.1°C) changes in temperature, induced by pulses 
from an infrared laser (see Methods). These stimuli are 
large enough to elicit responses but well below the 
threshold for pain avoidance . In Figure [S] we show 
the distribution pt(w) of phase velocities as a function 
of time relative to the thermal pulse. All of the worms 
were crawling forward at the moment of stimulation, so 
the initial phase velocities are distributed over a wide 
range of positive values. Within one second, the dis- 
tribution narrows dramatically, concentrating near zero 
phase velocity-the pause states described above. 

Arrival in the pause state is stereotyped both across 
trials and across worms. By analogy with conventional 
psychophysical methods [H , we can ask how reliably an 
observer could infer the presence of the heat pulse using 
the worm's response. We flnd that just measuring the 
phase velocity uj at single moment in time after the pulse 
is sufficient to provide ~ 75% correct detection of this 
small temperature change in single trials. 



COUPLING THE MODES AND STEERING THE 
WORM 

Our discussion thus far has separated the dynamics 
of the worm into two very different components: the 
{01,02} plane with its phase dynamics, responsible for 
crawling motions, and the mode 03 which is connected 
with large curvature turns. Because we found these 
modes as eigenvectors of a covariance matrix, we are 
guaranteed that the instantaneous amplitudes of the 
modes are not linearly correlated, but this does not mean 
that the dynamics of the different motions are completely 
uncoupled. We found the clearest indications of mode 
coupling between the phase in the {oi, 02} plane and the 
amplitude as at later times, which is illustrated by the 
correlation function in Fig[5|D. Intuitively, these correla- 
tions make sense, since the direction in which the worm 
forms its shape must depend on the phase of the wave 
passing along the worm's body at the moment that the 
large amplitude bending is initiated. This pattern of 
correlations is perturbed strongly by thermal stimuli. 

The fact that the correlations between phase and the 
turning mode are stimulus dependent implies that the re- 
sponse of the turning mode to thermal stimuli depends 
on the phase which the worm finds itself at the time of 
the stimulus. Perhaps some of the apparent random- 
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FIG. 6: Thermal responses, mode coupling and active steer- 
ing (see Methods for details), (a) The distribution of phase 
velocities pt{ijj) in response to a brief thermal stimulus. These 
data are taken from a collection of 13 worms, each stimu- 
lated with 20 repetitions of a AT = 0.1° C pulse. Note that 
within one second, the distribution becomes highly concen- 
trated near u) = Q, corresponding to the pause states identi- 
fied in Fig 5. In this experiment only, the worms were raised 
at a lower temperature leading to a lower average uj before 
the thermal stimulus, (b) Correlations between phase in the 
{01,02} plane and 03, C{t,t') = {5ai{t)5sm{4>{t'))). Shortly 
after the thermal impulse, the modes develop a strong anti- 
correlation which is distinct from normal crawling (c) Worm 
'steering'. A thermal impulse conditioned on the instan- 
taneous phase was delivered automatically and repeatedly, 
causing an orientation change €) in the worm's trajectory. In 
this example lasting 4 minutes, asynchronous impulses pro- 
duced a time-averaged orientation change (O) = O.Olrad/s 
(black), impulses at positive phase produced a trajectory 
with (O) = O.lOrad/s (blue), and impulses at negative phase 
produced (O) = — 0.12rad/s (red). This trajectory response 
is consistent with the early-time mode correlations seen in Fig 
6b. We found 13 out of 20 worms produced statistically differ- 
ent orientation changes under stimulated and non-simulated 
conditions while only 1 out of 20 worms responded in the 
same fashion when the phase was randomized. The asterisk 
marks the start of the trajectories. 



ness of turning responses thus is related to the fact that 
repeated thermal stimuli catch the worm at different ini- 
tial phases. As a preliminary test of this idea, we imple- 
mented our analysis online, allowing an estimate of the 
phase with a delay of less than 125 ms. We then deliver 
an infrared pulse when the phase falls within a phase 
window that corresponds to either dorsal- or ventral- 
directed head swings. To understand the predicted con- 
sequences of these pulses, consider an idealized example 
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in which the target phase is 7r/2; then the correlation 
structure in Fig[S}3 predicts a negative bias in 03 which 
should generate right hand turns. Similarly, if the tar- 
get phase is — 7r/2, then continued thermal stimulation 
should generate left hand turns. Note that the thermal 
pulse itself does not have a handedness, so that if the 
pulses are not synchronized to the state of the worm 
there should be no systematic preference for right vs. 
left handed turns. These predictions are confirmed in 
FigH;. 



DISCUSSION 

Our central result is that C. elegans motor behavior 
is vastly simpler than it could be, even when the or- 
ganism moves freely. Conceptually similar results have 
been obtained for aspects of motor control in humans 
and other primates, where postures or trajectories of 
limbs, hands or eyes are confined to spaces of low dimen- 
sionality despite the potential for more complex motions 
[20, 21, 22, 23J. For C. elegans itself, recent quantitative 
work has focused on simplifying behavior by matching 
to a discrete set of template behaviors, such as forward 
and backward motion of the center of mass , sinusoidal 
undulations of the body [H , or bends . Our results 
combine and generalize these ideas. Motor behaviors are 
described by projection of the body shape onto a small 
set of templates (the eigenworms), but the strengths of 
these projections vary continuously. The templates are 
sinuous, but not sinusoidal, because the fluctuations in 
posture are not homogeneous along the length of the 
worm. Our description of shape is intrinsic to the worm 
and invariant to the center of mass position, but mo- 
tion in shape space predicts the center of mass motion. 
There are discrete behavioral states, but these emerge 
as attractors of the underlying dynamics. Most impor- 
tantly, our choice of four eigenworms is driven not by 
hypotheses about the relevant components of behavior, 
but by the data itself. 

The construction of the eigenworms guarantees that 
the instantaneous amplitudes along the different dimen- 
sions of shape space are not correlated linearly, but the 
dynamics of the different amplitudes are nonlinear and 
coupled; what we think of as a single motor action always 
involves coordinating multiple degrees of freedom. Thus, 
forward and backward motion correspond to positive and 
negative phase velocity in Fig[3l but transitions between 
these behavioral states occur preferentially at particu- 
lar phases (data not shown). Similarly, H. turns involve 
a large amplitude excursion along 03, but motion along 
this mode is correlated with phase in the ({ai, 02}) plane, 
and this correlation itself has structure in time (FigOj). 
The problems of C. elegans motor control are simpler 
than for higher animals, but these nonlinear, coupled 



dynamics give a glimpse of the more general case. 

Perhaps because of the strong coupling between the 
turning mode 03 and the wriggling modes ai,a2, we have 
not found an equation of motion for 03 alone which would 
be analogous to Eq's (|4]) for the phase. Further work is 
required to construct a fully three dimensional dynamics 
which could predict the more complex correlations such 
as those in Fig [61d. Turning should emerge from these 
equations not as another attractor, but as an 'excitable' 
orbit analogous to the action potential in the Hodgkin- 
Huxley equations or to recent ideas about transient dif- 
ferentiation in genetic circuits jlJl. A major challenge 
would be to show that the stochastic dynamics of these 
equations can generate longer sequences of stereotyped 
events, such as pirouettes [25]. 

We have shown that a meaningful set of behavioral 
coordinates can uncover deterministic responses. A re- 
sponse might seem stochastic or noisy because it depends 
on one or more behavioral variables that are not being 
considered. In our experiments, nonlinear correlations 
among the behavioral variables suggest that some of the 
randomness in behavioral responses could be removed if 
sensory stimuli are delivered only when the worm is at 
a well defined initial state, and we confirmed this pre- 
diction by showing that phase-aligned thermal stimuli 
can 'steer' the worm into trajectories with a definite chi- 
rality. A crucial aspect of these experiments is that the 
stimulus is scalar — a temperature change in time has 
no spatial direction or handedness — but the response, 
by virtue of the correlation between stimulus and body 
shape, does have a definite spatial structure. The align- 
ment of thermal stimuli with the phase of the worm's 
movement in these experiments mimics the correlation 
between body shape and sensory input that occurs as 
the worm crawls in a thermal gradient, so the enhanced 
determinism of responses under these conditions may be 
connected to the computations whic h g enerate nearly 
deterministic isothermal tracking 1^, 26| . 

More generally, all behavioral responses have some 
mixture of deterministic and stochastic components. In 
humans and other primates, it seems straightforward 
to create conditions that result in highly reproducible, 
stereotyped behaviors, such as reaching movements 27 1. 
At the opposite extreme, bacterial motility is modulated 
in response to sensory inputs, but these responses seem 
fundamentally probabilistic [28]. Some of these differ- 
ences may result from the physical nature of sensory 
stimuli in organisms of vastly different size 2^ 3^ , but 



some of the differences may also result from differences 
of strategy or available computational power. The more 
stochastic the response, the more challenging it is to 
characterize behavior quantitatively and to link behav- 
ior with underlying molecular and neural components, 
as is clear from recent work on Drosophila olfaction (see, 
for example, Ref [sij). We hope that our approach to 
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the analysis of behavior may help to uncover more de- 
terministic components of the sensory-motor responses 
in other model organisms. 

More than forty years of work on C. elegans has led 
to a fully sequenced genome [12] and to the complete 
wiring diagram of the nervous system (33| . Significant 
steps have been made toward the original dream [3411 
of connecting genes, neurons, and behavior [sl, l36l |. 
Nonetheless, with the advances in molecular, cellular, 
and circuit analyses, our ability to probe the mechanisms 
which generate behavior substantially exceeds our ability 
to characterize the behavior itself. Perhaps our work 
provides a step toward addressing this imbalance. 



METHODS 

Tracking microscopy. The imaging system con- 
sists of a Easier firewire CMOS camera (A601f, Easier, 
Ahrensburg, Germany) with 4x lens (55-901, Ed- 
mund Optics, Earrington, NJ) and a fiber optic trans- 
illuminator (DC-950, Dolan-Jenner, Eoxborough, MA) 
mounted to an optical rail (Thorlabs, Newton, NJ). 
The rail is attached to a XY translation stage (Deltron, 
Bethel, CT) which is driven by stepper motors (US Digi- 
tal, Vancouver, Washington). The stage driver is a home 
made unit utilizing a SimpleStep board (SimpleStep, 
Newton, NJ) and Gecko stepper motor drivers (Geck- 
odrive, Santa Ana, CA). Image acquisition, processing, 
and stage driver control was done using Lab VIEW (Na- 
tional Instruments, Austin, TX). Images of worms were 
isolated and identified using the image particle filter. A 
raw unprocessed JPEG image and a filtered process bi- 
nary PNG image were written to the hard drive at rates 
up to 32 Hz. Concurrently at 4 Hz, the center of mass 
of the worm was calculated and the distance from the 
center of the field of view in pixels was computed. An 
error signal was then calculated via a coordinate trans- 
formation between the camera reference frame and the 
translational stage reference frame and the XY stage was 
moved to center the worm in the field of view. 

Worm preparation. The C. elegans strain, N2, 
was grown at 20° C and maintained under standard con- 
ditions Before each experiment, excess moisture 
from NGM assay plates (1.7% Eacto Agar, 0.25% Eacto- 
Peptone, 0.3% NaCl, 1 mM CaCls, 1 mM MgS04, 25 
mM potassium phosphate buffer, 5 /ig/mL cholesterol) 
was removed by leaving them partially uncovered for 1 
hr. A copper ring (5.1-cm inner diameter) pressed into 
the agar surface prevented worms from crawling to the 
side of the plate. Young adults were rinsed of E. coli by 
transferring them with a worm pick from OP50 bacterial 
food plates into NGM buffer (same inorganic ion concen- 
tration as NGM assay plates) and letting them swim for 
1 minute. Worms were transferred from the NGM buffer 



to the center of the assay plate (9-cm Petri dish). The 
plates were covered and tracking began after 1 minute 
and lasted no longer than 60 minutes. In the rare cases 
where worms stopped moving before the completion of 
the run, the data were excluded. 

Eigenviforms. Images of worms captured by the 
worm tracker were processed using MATLAE (Math- 
works, Natick, MA). Cases of self- intersection were ex- 
cluded from processing. Images of worms were thinned 
to a single-pixel-thick backbone. A spline was fit through 
these points and then discretized into 101 segments, 
evenly spaced in units of the backbone arclength. The 
N = 100 angles between these segments were calcu- 
lated and an overall rotation mode was removed by 
subtracting X] ^('^(O)/-^ from each angle. The shape 
covariance matrix C(s,s') = {{9{s) - (9)) (6'(s') - (9))) 
was constructed from 9 freely crawling worms, sam- 
pled at 4 Hz, for a period of 30 minutes. Each eigen- 
worm u^(s) is an eigenvector of the covariance matrix 
(^(s, s')'"m('5') = ^fj-U^is). The fractional variance 
captured by K eigenvectors is thus aj^ — J2^=i ^fj./^'^y 
where — ^^'^ total variance of the measure- 

ments. The same eigenworms shown in Fig. [5] were used 
throughout the various analysis reported in the paper. 
The worm's phase was defined as (j) ~ tan^^ (— a2/ai) 
where ai and 02 were both normalized to unit variance. 

Equations of motion. For the analysis of phase dy- 
namics we sampled the worm shape at 32 Hz. Data for 
the construction of the equations of motion came from 
12 worms, 5 trials per worm, with 4000 frames per trial. 
We also filtered each mode time series through a low-pass 
polynomial filter so that for each frame (26 < m < 3974), 

= J2lt-.25Yl'j=oPo(^ ~ ^^^^^ fe} ^"^^ the 
best-fit polynomial coefficients. Mode time derivatives 

were calculated using derivatives of the polynomial filter. 
None of our results depend critically on the properties of 
the filter. The Langevin equations governing the phase 
dynamics are shown Eq. ([4]) and we learn the functions 
{F{(f>,uj),a{(j),uj)} directly from the time series H, 3^. 
By construction {a[(j){t),(jj{t)]r]{t)) — and therefore 
the optimal rms estimate of F{(j), lu) is the conditional 
mean {Lu\uj,(j)). We estimate F by assuming a functional 
expansion F{uj,(j)) = Z]f„=-5 Z]p=o '^m^'^e"""'^, where 
the model parameters {aj^} were determined by mini- 
mizing the rms error = J2t ('^(*) ~ ^['^(*); 
training data (90%) and the hyperparameters {mmax = 
= 5} were chosen to minimize error on held- 
out data (10%). Once F is known we can determine 
the noise in the system; we normalize (77^) so that 
a^{uj,(f>) — {{lu ~ F{uj,(f>)) \uj,(f>). The attractors con- 
tained within our derived dynamics were obtained by 
evolving initial conditions spanning the sampled {to, (j)} 
plane for long times (93.75s ^ 47 cycles). In the deter- 
ministic dynamics all trajectories evolve to one of four 
asymptotic states and we observed no switching. 
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Thermal impulse response (experiment). 

Worms were prepared as described earlier. A collimated 
beam with a 1/e diameter of 5.6 mm (standard stim- 
ulus) or 1.5 mm (painful) from a 1440 nm diode laser 
(FOL1404QQM, Fitel, Peachtree City, GA) was posi- 
tioned to heat the area covering the worm. The diode 
laser was driven with a commercial power supply and 
controller (Thorlabs, Newton, NJ). Power and duration 
of the beam was controlled through software using 
Lab VIEW. For each worm, 1000 seconds of data was 
collected in cycles of 50 seconds. 12.5 seconds into each 
cycle the laser was turned on for a duration of 75 ms 
at 150 mW (standard) or 250 ms at lOOmW (painful). 
The temperature increase caused by the laser pulses 
was measured using a 0.075mm T-type thermocouple 
(coco-003. Omega, Stamford, CT) placed on the surface 
of the agar and sampled with a thermocouple data 
acquisition device (USB-9211, National Instruments). 
For each measurement, 60 trials of 30 s cycles were 
averaged. The temperature increase was calculated 
by subtracting the maximum temperature (recorded 
immediately after the laser pulse) from the baseline 
temperature (recorded 9 s after the laser pulse). The 
temperature increase for the standard pulse was 0.12°C 
and the increase of the painful pulse was 0.73°C. 

Thermal impulse response (analysis). In Fig 

6a, the time-dependent probability density pt(a;) was 
smoothed before the onset of the impulse with a gaus- 
sian low-pass filter of size 0.19 s in the t direction and 
0.17cycles/s in the u direction. In Fig 6b the correlation 
function C(t,t') = ((sin^(t) - (sin./))) (03(^0 - (03))) 
was calculated as follows. Far from the time of the 
impulse (frames 800 to 1574, impulse on frame 400), 
we expect time-translation invariance Cpost{t,t') = 
9{t - t') + €post(i,i') where 5((A) = (C(i, j))i-j=A 
is the true correlation function and ^post(^, i') charac- 
terizes statistical error. Similarly in a time window 
aroimd the impulse (frames 24 to 800), Cstim{t,t') = 
g{t — t') + ^stim(i,i')- However, the thermal im- 
pulse breaks this invariance and ^siim(i.i') contains 
both sampling fluctuations and stimulus-dependent cor- 
relation dynamics. To separate these eff'ects we 
use singiilar value decomposition to compare ^post 
and ^stim- We write each matrix Cpost/stim(i, *') = 

Tilt" ^post/stim(^i ^")'5'post/stim(^"j '^")^ost/stim('^"j ^ud 

find that only two singular values of ^stim arc signifi- 
cantly larger then ^post- We then reconstruct the two- 
point function around the stimulus as Cstim{t,t') = 

9{t - 1') + Et"=i u,tUt, i")5stim(t", t")v,,ut", t'). 

Thermal steering. Preparation of worms and in- 
strumentation were the same as described for the ther- 
mal impulse response. However, instead of process- 
ing worm images off-line, real-time calculation of the 
eigenworms and shape phase </> was done using custom 
dynamic-linked image processing libraries written in C 



along with supporting Lab VIEW code. The modes were 
computed as previously described except that the spline 
interpolation algorithm was replaced with a Hermitian 
interpolation algorithm to reduce the processing time. 
The processing time was short enough to simultaneously 
track and calculate modes at 8 Hz. For phase dependent 
measurements, the laser was fired when the worm was 
moving forward and (p fell within a prescribed interval 
(width 1 radian). The laser pulse (150 mW) lasted for 
75 ms and caused a temperature increase of 0.12° C. For 
each run a pair of triggering phase windows (0 to -1, and 
2.1 and 3.1 radians) corresponding to the dorsal- and 
ventral-directed head swing was used. The sequence of 
each run started with a 5 minute period of no stimulus 
followed by the pair of phase dependent stimuli. The or- 
der of each pair of stimulus conditions was switched for 
each successive run. For the randomized pulse control 
experiments, the laser was fired with a uniform phase 
probability, but with conditions that restricted the fir- 
ing interval to be longer than 2 seconds. 

Steering and turn identification. The time- 
average change in orientation of the worms path, (0) 
(rad/s), was calculated from the angular changes be- 
tween the positions of the center of mass of the worm 
during forward runs of at least 4 s in length. Given posi- 
tions (ri, r2, rs, . . . , rjv), the angles between connecting 
segments (r2 — n, rs — r2, . . . , rjv — r(]v-i)) were calcu- 
lated. (0) was calculated in intervals of 10 s. Since the 
distributions were Gaussian (data not shown) with sim- 
ilar variance, we used the Student's t-test to determine 
if the values of & under thermal stimulation were sig- 
nificantly different than the control {p < 0.05). Since 
we were interested in the change in orientation dur- 
ing forward motion we excluded trajectory data that 
contained large turns or reversals along with angular 
changes greater than 7r/4 radians. These events were 
automatically detected by measuring the compactness of 
the worm shape. Compactness was calculated by mea- 
suring the longest distance between two points in the 
worm shape (also known as the max feret distance) and 
normalizing this with the maximum value for the entire 
data run. Turns were flagged when the compactness fell 
below 0.6. 
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